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ABSTRACT 

We study the structure and evolution of “quasi-stars’, accreting black holes embedded within 
massive hydrostatic gaseous envelopes. These configurations may model the early growth of 
supermassive black hole seeds. The accretion rate on to the black hole adjusts so that the 
luminosity carried by the convective envelope equals the Eddington limit for the total mass, 
M» + Mgy ~ M:. This greatly exceeds the Eddington limit for the black hole mass alone, 
leading to rapid growth of the black hole. We use analytic models and numerical stellar structure 
calculations to study the structure and evolution of quasi-stars. We show that the photospheric 
temperature of the envelope scales as Tp, « Moy °°, and decreases with time while the 
black hole mass increases. Once Tph < 10* K, the photospheric opacity drops precipitously and 
Tph hits a limiting value, analogous to the Hayashi track for red giants and protostars, below 
which no hydrostatic solution for the convective envelope exists. For metal-free (Population III) 
opacities, this limiting temperature is approximately 4000 K. After a quasi-star reaches this 
limiting temperature, it is rapidly dispersed by radiation pressure. We find that black hole 
seeds with masses between 10° and 104 Mo could form via this mechanism in less than a 
few Myr. 
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1 INTRODUCTION 


The formation mechanism for supermassive black holes in galactic 
nuclei remains unknown. Variations on most of the formation chan- 
nels identified by Begelman & Rees (1978) — which include insta- 
bilities in clusters of stars or stellar remnants and the collapse of 
supermassive stars or massive discs — are still under consideration 
today (Umemura, Loeb & Turner 1993; Shibata & Shapiro 2002; 
Freitag, Giirkan & Rasio 2006; Lodato & Natarajan 2006). A more 
recent idea holds that supermassive black holes result from sustained 
accretion on to, or mergers of, the remnants of Population III (Pop 
IID (i.e. metal-free) stars (Volonteri, Haardt & Madau 2003), at least 
some of which seems likely to be massive, short-lived progenitors of 
stellar mass black holes (Carr, Bond & Arnett 1984; Bromm, Coppi 
& Larson 1999; Abel, Bryan & Norman 2002; Heger et al. 2003; 
Tumlinson, Venkatesan & Shull 2004). 

No current observation directly constrains the different classes of 
model for black hole formation. However, there are clues. The ex- 
istence of massive black holes in quasars at redshifts z > 6 strongly 
suggests that these black holes, at least, started forming at redshifts 
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high enough (z ~ 20) to predate extensive metal enrichment by 
the first stars. Motivated by this, we develop in this paper a model 
proposed by Begelman, Volonteri & Rees (2006b hereafter BVR), 
in which supermassive black holes form, not from Pop III stars 
themselves, but rather from the evolution of a new class of Pop III 
objects that might form in metal-free haloes too massive to yield 
individual stars. In outline, the BVR model envisages a three-stage 
process for black hole formation. First, gas in metal-free haloes 
with a virial temperature above T ~ 10* K flows towards the cen- 
tre of the potential as a result of gravitational instabilities, form- 
ing a massive, pressure-supported central object. Nuclear reactions 
may start, but the very high infall rate continues to compress and 
heat the core, precluding formation of an ordinary star. Eventually, 
when the core temperature attains T ~ 5 x 108 K, neutrino losses 
result in a catastrophic collapse of the core to a black hole. We 
dub the resulting structure — comprising an initially low-mass black 
hole embedded within a massive, radiation-pressure-supported en- 
velope — a quasi-star. Initially, the black hole is much less massive 
than the envelope. Over time, the black hole grows at the expense 
of the envelope, until finally the growing luminosity succeeds in 
unbinding the envelope and the seed black hole is unveiled. The 
key feature of this scenario is that while the black hole is em- 
bedded within the envelope, its growth is limited by the Edding- 
ton limit for the whole quasi-star, rather than that appropriate for 
the black hole mass itself. Very rapid growth can then occur at 
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Figure 1. Schematic illustration of the quasi-star structure that we consider 
in this paper. A seed black hole of mass Mgy accretes gas from a massive, 
radiation-pressure-supported envelope at a rate set by the conditions outside 
the Bondi radius. The luminosity liberated by the accretion process is trans- 
ported convectively in the inner regions of the envelope, with a transition 
to a radiative zone once convection becomes inefficient. In this paper, we 
consider isolated, spherically symmetric models of quasi-stars, but the more 
physical situation would also include partial rotational support leading to 
flattening of the quasi-star, and on-going disc accretion at a fraction of a 
Solar mass per year. 


early times, when the envelope mass greatly exceeds the black hole 
mass. 

The structure of quasi-stars, illustrated in Fig. 1, resembles that 
of more familiar objects. The outer regions have qualitative simi- 
larities to Thorne & Zytkow (1977) objects, although in our case 
the luminosity derives entirely from black hole accretion. The ac- 
cretion physics is related to the ‘hypercritical’ regime of accretion 
(Begelman 1979; Blondin 1986), which may be observed in X-ray 
binaries such as S8433 (Begelman, King & Pringle 2006a). How- 
ever, these are very rough analogues. The large difference between 
Pop III opacity and that of present day, metal-enriched gas, to- 
gether with the dominant role of radiation pressure and the presence 
of on-going accretion in quasi-stars, results in a unique quasi-star 
structure. 

In this paper, we study the structure and evolution of quasi-stars. 
Using Pop III opacities, we construct simplified analytic models 
which we compare with numerical integrations of the stellar struc- 
ture equations. Our model quasi-stars are assumed to be spherically 
symmetric and in hydrostatic equilibrium. Real quasi-stars, if they 
exist, may not obey these strictures, and thus our results cannot be 
definitive. Rather, our calculations are intended to provide a first 
estimate of the maximum black hole mass that this mechanism can 
yield, together with a determination of the photospheric tempera- 
ture and luminosity of quasi-stars that can be used to assess their 
detectability with next-generation observatories such as the James 
Webb Space Telescope. 


2 ANALYTIC CONSIDERATIONS 


We consider a black hole of mass Mgy embedded within an envelope 
of mass M, >> Mpy. The luminosity of the quasi-star is generated 
exclusively by black hole accretion, whose rate depends upon the 
central conditions in the envelope. To proceed analytically, we first 
develop simple scaling relations for the envelope structure and black 
hole accretion rate. 


2.1 Envelope structure 


Quasi-star envelopes with masses > 10° Mo are supported primarily 
by radiation pressure. Since the luminosity carried by the envelope 
must equal the Eddington limit for the total mass, the flux carried 
by radiative diffusion in the envelope’s interior is only a fraction 
M(<r)/M, of the total, where M(< r) is the mass enclosed within r. 


Thus, quasi-star envelopes are strongly convective, and their struc- 
tures resemble n = 3 (y = 4/3) polytropes. The most accurate 
approach would be to model the envelope as a ‘loaded polytrope’ 
(Huntley & Saslaw 1975), with the black hole treated as a central 
point mass, but for M, >> Mgy the standard Lane—Emden solutions 
suffice. Defining m, = M,/Mọ , the ratio of gas to radiation pressure 
is uniform throughout the convective zone. We find 


2 =7.1mz'”, (1) 
Pr 

where we have assumed a mean mass per particle u ~ 0.6m, valid 
for fully ionized regions. 


2.2 Energy source 


The central regions of n = 3 polytropes have approximately uniform 
densities (o.), temperatures (T<) and pressures (pe = aT$ /3). The 
boundary conditions for black hole accretion are therefore similar 
to those of Bondi (1952) accretion, or more specifically the gen- 
eralization for optical thick flows given by Flammang (1982). The 
adiabatic accretion rate is 
4r (GMpn} Pe 
TE 
where ce = (4p /3p.)!/? is the adiabatic sound speed. Bondi’s solu- 
tion assumes that all gravitational binding energy liberated during 
accretion is advected into the hole, but this is unrealistic in the pres- 
ence of even a small amount of rotation. Provided that the specific 
angular momentum at the Bondi radius fow exceeds the specific 
angular momentum lms of the marginally stable circular orbit, we 
expect the flow to be rotationally supported. In this case, a geomet- 
rically thick accretion disc will form around the hole, within which 
angular momentum transport is required in order for accretion to 
occur. Although the efficiency of such a disc (and how the energy 
outputis partitioned between radiation and mechanical work) is hard 
to calculate, it is reasonable to assume that it will depend primarily 
on the depth of the black hole potential. We write the luminosity 
as Ley = €Mpyc?, where € ~ O(0.1) is the efficiency of energy 
output! and Mpy is the actual accretion rate. For simplicity, we take 
€ to be a constant. In the absence of an efficient exhaust, such as 
a jet or evacuated funnel, this energy must be carried beyond the 
Bondi radius convectively, by the accreting gas itself. Given that 
the convective flux density may not exceed ~ p.c., we conclude 
that the accretion rate will be reduced below the Bondi value by a 
factor of the order of €~!(c./c)? « 1 (Gruzinov 1998; Blandford 
& Begelman 1999; Narayan, Igumenshchev & Abramowicz 2000; 
Quataert & Gruzinov 2000). Unless € « 0.1 — for example, in the 
case where Igy < lms — accretion is suppressed by a large factor, 
of the order of 10* for the models we consider. Combining these 
results we adopt the expression 


Mbondi = 


, (2) 


Len = 4nG’a Mgao py” (3) 


for the accretion luminosity. The parameter œ < 1 accounts for en- 
ergy sinks within the Bondi radius: inefficient convection, presence 
of outflows, etc.; as well as any inefficiency of angular momentum 
transport. Small values of œ imply a reduced energy supply to the 


1 We note that € is expected to vary by factors of the order of unity with the 
black hole spin. In our models, however, we are envisaging a black hole that 
grows in mass by several orders of magnitude. If the accretion flow retains a 
similar character during this growth, we expect that the spin of the hole will 
rapidly attain a limiting value, after which further changes will be small. 
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quasi-star. The standard result for an n = 3 polytrope (e.g. Hoyle & 
Fowler 1963), 


pe = 1.3 x 10-4mz!?T3 gem, (4) 


where Te = T./10° K, coupled with the equation of state for radia- 
tion in local thermodynamic equilibrium (LTE), allows us to express 
Legu in terms of the central temperature: 

Leu = 6.6 x 10am ymn Te? ergs |, (5) 


where MBH = Mpu/Mo. 


2.3 Scaling laws 


Suppose that the convective zone encompasses nearly the entire 
mass and radius of the envelope, allowing us to treat M, and R, as 
constant within the radiative layer. We will check the validity of this 
assumption in Section 2.4.2. We expect Lgu to be very close to the 
Eddington limit at the transition between the convective zone and 
the outer layer where radiative diffusion carries all the flux, and set 
Lg = 4n GM,,c/k 1, Where Ky is the opacity at the transition radius. 
As we will see below, ky is close to the electron scattering opacity, 
Kes = 0.35 cm? ots hence Š = Kir/Kes ~ O(1). We therefore write: 


Lpy = 1.4 x 10° 2,,@7'm, erg gt, (6) 


where £y œ 1 is the Eddington factor at the transition. Equating the 
two expressions for Lgy, we obtain 


Ty = 1.4 x 10 am mK. (7) 


We are interested in systems with m, ~ O(100 mgn), so T, typi- 
cally lies in the range 10° to 10° K and electron scattering opacity 
dominates in the quasi-star interior. At these temperatures, the rates 
of energetically significant thermonuclear reactions are negligible, 
and can be safely ignored”. We note that this differs from the case of 
envelope accretion on to neutron stars, where the presence of a hard 
surface results in high temperatures and potentially significant nu- 
clear reactions within the hydrostatic region of the flow. Although 
we do not model the non-hydrostatic region of quasi-stars in any 
detail, in our case high temperatures are attained in the immediate 
vicinity of the black hole. However, even in this region the neglect 
of nuclear reactions is justified, first because black hole accretion is 
energetically much more efficient than fusion, and secondly because 
on-going accretion limits the time-scale over which inflowing gas 
is exposed to high T. 

The polytropic relations also allow us to estimate the radius of 
the quasi-star: 


R, = 5.8 x 10°m!’ Tg ' cm, (8) 


(Hoyle & Fowler 1963) and therefore its photospheric (effective) 
temperature, 


Ton = 8.5 x 10L RTI’ K. (9) 
Inserting our estimate for T., we obtain 

R, = 4.3 x 1003 nRa mgm, cm (10) 
and 

Tyr = 1.0 x 1004/0 R Vym im K, (11) 


2 In the most massive quasi-stars, the central temperature may be high enough 
(a few million K) to initiate lithium burning. This is energetically negligible, 
and although the presence or absence of lithium does affect the opacity, the 
effect is small for the photospheric temperatures and densities of interest 
here. 
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Thus, quasi-stars should have radii of the order of 107—10° au and 
temperatures of a few thousand degrees. 


2.4 Radiative layer and photosphere 


In the outer layers of the quasi-star, convection is unable to trans- 
port the total luminosity and radiative diffusion takes over as the 
dominant energy transport mechanism. To estimate the transition 
temperature between the convective and radiative zones, we note 
that the maximum flux that can be transported convectively is 


Foon, max = Bpc,, (12) 


where c, = (p/p)! is the local (isothermal) sound speed and B < 
1 is an efficiency factor. Equating this to the photospheric flux 
acT on, /4 (again assuming a narrow radiative zone), we obtain 
3 c 

4 4 

p= ap cm (13) 
This result is equivalent to the condition that radiation be trapped 
in the convective cells as they rise. Since B < 1 and c/c, > 1, we 
deduce that Ty >> Tpn. Expressing p and p in terms of Ty (using the 
polytropic relations and the equation of state), we obtain 


Ta = 3186Om ha. (14) 


Since we will later show that T,, cannot drop below a few thousand 
degrees, we conclude that Ty is well above 10* K. The import of 
this is that at the densities and temperatures likely to apply near 
the transition (for m, > 10+), bound—free opacity can be important 
but will not elevate the Rosseland mean opacity above Thomson 
scattering by a large factor. 

In order to determine the structure of the radiative layer, we need 
to know the behaviour of the opacity. For the purpose of our analytic 
estimates, we adopt a simple phenomenological form for the opac- 
ity, based on the Pop III opacity tables of Mayer & Duschl (2005, 
hereafter MD05). We note that, for temperatures below 10* K and 
densities <10~° g cm~?, the Rosseland mean opacity depends much 
more sensitively on temperature than on density. This can be seen 
clearly in fig. 4 of MDOS. We therefore write 


Ko 
K(T) = LET (15) 
This expression does not capture the possible contribution of bound— 
free opacity above 10* K, but it does mimic the extremely steep 
decline in opacity towards lower temperatures. A crude fit gives 
s Z 13, To © 8000 K and ko & Kes. 

The advantage of using an opacity that is solely a function of 
temperature is that one can combine the equation of hydrostatic 
equilibrium with the radiation diffusion equation and integrate to 
obtain the total (radiation + gas) pressure as a function of tempera- 
ture (see e.g. Cox 1968, chapter 20, or any stellar structure textbook 
for a discussion of this technique). Defining the Eddington factor 
associated with opacity Ko as 


Lguko 
bo = ———_ ., 16 
0 = 4nGM,c (o) 
we obtain 
da |T- T4 R(T -T 
p(T) = 30 7 at ( z 7 + Pph, (17) 
0 SS 


where ppn is the total pressure at the photosphere. We use a standard 
model for the atmospheric structure, in which the radiation pressure 
depends on optical depth as p,(t) = (aT}, /6) (1 + 31/2) and the 
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temperature is constant between the true (t = 0) surface and the 
photosphere at t = 2/3 (see e.g. Mihalas 1978 or most stellar struc- 
ture textbooks). Integrating the equation of hydrostatic equilibrium 
and noting that the surface gravity can be written as g = ko aT Sy / 
(4 £o), we have 


2 g aTi a TNS 
= F = 1+8 +| = ; 18 
Poh = 3 Ta) 6 Os ae Ton k 
The gas pressure is then given by 
(T) = aT* 
Pg =P 3 
4a T* Ti 
= 1-¢ r 
300 { o ( 4 8 
Ty s+4 
= RY a ( 8 Lie ry] ; (19) 


for temperatures Tpn < T < Ty. 

The ratio of gas pressure to radiation pressure must be continuous 
across the boundary between the convective zone and the radiative 
layer. The matching condition is, therefore, 


3 p(T) 
ai, 


=7.1mz'/?. (20) 


To simplify matters, we note that, since Te >> Tp, and s > 4, we 
can neglect the terms in p(T) proportional to Toh and T**. Solving 
the matching condition for £9, we obtain 


1+ (s+ 4)/[2s — 4)] [To/(Ta° (Tir/Ton) 
bo = = i 21) 
1 +7.1m, 


2.4.1 Opacity crisis 


A necessary condition for a model to be physically realistic is that 
the flux not exceed the Eddington limit at the transition radius, 
fy < 1. At the transition radius, radiation is effectively trapped, 
and any super-Eddington flux would be efficiently converted into 
bulk kinetic energy. Hence, if this condition is violated we expect 
the outer layers to be blown off by radiation pressure, ultimately 
dispersing the quasi-star. For our approximate opacity, which is a 
monotonically increasing function of T, this translates to £9 < 1 + 
(To/T.)’. Since the convective envelope is dominated by radiation 
pressure, we can write this condition as 


s s—4 s 
s +4 T Ty T 
DAE 0 : 2Fim +2). (22) 
2(s 4) Tr Tn Tr 


Noting that the second term on the right-hand side of this equation is 
negligible, we substitute for T, using equation (14) to obtain a lower 
limit on Tp, 


s+4 
PP > 7.7 x 10 (= iz +) pP mST. (23) 


For s > 4, this lower limit is very insensitive to the convective 
efficiency £ and the envelope mass M,, and is roughly proportional 
to To. For example, adopting s = 13 and normalizing 6,m, and To 
to 0.1, 10* and 8000 K, respectively, we have 


Ty, > 4500 pois, aay pe K. (24) 
The corresponding lower limit on Ty shows similar behaviour: 
Ta > 55000 Bo 18/113 ra eas K. (25) 


The first of these lower limits is the most important result of 
the analytic part of this paper. It represents a floor to the photo- 


spheric temperature of quasi-stars, analogous to the ‘Hayashi track’ 
(Hayashi 1961; Hayashi & Hoshi 1961), which limits the temper- 
atures of red giants and convective protostars. Our analysis differs 
from Hayashi’s model in that our convective envelopes are radiation- 
pressure dominated and therefore resemble n = 3 polytropes, rather 
than the n = 3/2 polytropes appropriate to gas-pressure-dominated 
convection. 


2.4.2 Validity of approximations 


The analytic model described above assumes that the growth of the 
black hole within the envelope can be described by a sequence of 
hydrostatic solutions. This quasi-static approximation is justified by 
the ordering of time-scales. Adopting representative values of the 
mass, radius and radiative efficiency (m, = 10°, mpy = 10°, R, = 
105 cm, € = 0.1), the time-scale on which hydrostatic equilibrium 
is established, 


GM, 
~ 108s (26) 


tayn me 
* 


is much shorter than the time-scale on which the black hole grows: 


M 
teow = = ~ 10” s, (27) 
Mren 
Similarly, the time-scale for attaining radiative equilibrium 
TAR, 


c 


frad ~ ~ 1078, (28) 


where t is the optical depth through the radiative layer of width 
AR,,, is enormously shorter than any of the evolutionary time-scales 
of the system. 

In addition to the quasi-static assumption, we have taken the ge- 
ometric thickness and mass of the radiative layer to be negligible 
compared to R, and M,, respectively. How good are these assump- 
tions? 

First, consider the mass of the radiative layer. If the layer is ge- 
ometrically thin and the gravity g is approximately constant across 
it, then we have 


AM, _ 4nR*aT? 

S ae (29) 
M, 3GM2 

from the equation of hydrostatic equilibrium. Using equations (10), 

(11) and (14) to eliminate R,, Ty and mpy in favour of Tpn, we obtain 

(for £a = K = 1) 


40/9 
AM, -8/9 2/9 Ton 
zx 0.2 ——_ g 30 
“M, Boy Mas 4500K GO) 


where 6_; = 6/0.1. Equation (30) indicates that the assumption of 
constant mass in the radiative zone is only marginally self-consistent 
when T,y is close to its minimum value. The approximation improves 
at larger Tp. 

The geometric thickness of the radiative layer is dominated by the 
region close to the transition temperature. The opacity is therefore 
close to Ko ~ Kes and we may approximate 


AR,  167R,acT; 
R, BKes Ptr 


T —10/9 
= 0.78m ae ph 3] 
Boy Ma 4500 K GN 


Evidently, assuming geometrical thinness in the radiative layer is an 
even poorer approximation than assuming constant enclosed mass. 
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In fact, our numerical models show that even the scaling in equa- 
tion (31) — obtained assuming AR,/R, « 1 — is not valid. The 
radiative layer thickness generally decreases with T n- 


2.5 Stability 


The interior conditions of quasi-stars are hot enough that the opacity 
is dominated by electron scattering, yet too cool for nuclear reactions 
to occur. Accordingly, we do not expect any of the stellar instabilities 
that depend upon complex opacities or nuclear reaction rates to 
afflict the interior of quasi-stars (we note later the possibility of 
an instability near the surface due to the presence of locally super- 
Eddington zones). A more serious concern — given that quasi-stars 
are highly radiation-pressure dominated — is dynamically instability. 
In anon-rotating model, relativistic effects raise the critical y, below 
which instability occurs, to 


4 
Verit = 3 + Ô Verits (32) 
where ôy cit ~ GM/Rc? (e.g. Shapiro & Teukolsky 1983). This must 
be compared to the actual y in the radiation-dominated interior of 
the quasi-star. For p,/p < 1, the first adiabatic exponent is given by 
4 4 1p, 


MEg Pyes tae (33) 
If the pressure-weighted integral over the star of dy < dy cit, insta- 
bility is possible. 

The region of the quasi-star interior to the Bondi radius is not 
in hydrostatic equilibrium, so it is meaningless to apply the above 
condition there. For an estimate, we evaluate the above expressions 
at the Bondi radius using the expressions for the interior structure 
of the quasi-star. We find that the change in the stability boundary 
is 


2 
fig (2) = 4 x 10° mma, (34) 


while the actual y exceeds 4/3 by an amount, 
ôy > 1.2m". (35) 


We conclude that instability is possible if m, is large enough. For 
a = 0.01 and mgu = 10°, for example, instability is possible for 
m, > 3 x 106. 

We caution that this analysis is far too simple to describe the 
actual situation. Globally, a robust expectation is that quasi-stars 
will be rapidly rotating as a consequence of their formation from 
rotationally supported gas (BVR). Unfortunately, this does not help 
us in estimating the rotation rate at the Bondi radius, as the latter 
depends upon the very uncertain role of convection in redistributing 
angular momentum. However, any significant rotation would help 
stabilize the structure against dynamical instability. 


3 NUMERICAL MODELS 


Although our analytic models ought to yield a reasonably accu- 
rate picture of the generic features of quasi-star envelopes, they are 
only marginally self-consistent for photospheric temperatures ap- 
proaching the floor value. Moreover, real Pop III opacities, which 
are functions of density as well as of temperature, have considerably 
more structure than the analytic fit employed in the analytic models, 
and require a numerical treatment. 
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3.1 Equations 
We model quasi-stars as static, spherically symmetric objects in 
thermal equilibrium. The equation of hydrostatic equilibrium is 
dp GM(r) 
dr r? 


where the mass enclosed within radius r is given by 


p, (36) 


dM 5 
a (37) 
dr 


At the Bondi radius, which for our numerical integrations we define 
as 


Rona = PREP, (38) 
De 

we set M(Rponai) = Mpu. The gas mass within the Bondi radius is 

negligibly small compared to the black hole mass, so ignoring the 

non-hydrostatic region in our structure calculations should be an 

excellent approximation. The equation of state is a mix of gas and 

radiation pressure: 


P = Pe + pr = T +-aT*, (39) 


where k is the Boltzmann constant and u is the mean molecular 
weight. We set u = 0.6m,, appropriate for a fully ionized gas, and 
ignore the variation of u due to partial ionization near the quasi-star 
photosphere. This is a very good approximation, since the neutral 
hydrogen fraction generally does not exceed a few per cent. 

The temperature gradient is determined by the relative amount 
of energy carried by convection (Fon) and radiation (Fraa) at each 
radius. Since quasi-stars are powered solely by black hole accre- 
tion, the luminosity Ley = 47tr?(Feon + Fraa) is constant with ra- 
dius. To determine the boundary of the convective zone, we use the 
Schwarzschild criterion for convective stability,’ 


dT rad dT a 
dr dr 
where the radiative gradient dT ,aq/dr is the gradient of temperature 


if all the energy is transported by radiative diffusion: 


> 0, (40) 


dTraa 3 Lgu pK 
= ; 41 
dr 4ac 4nr? T? 4D 
and the adiabatic gradient: 
dTa r2—-1\Td 
t= (2 2, (42) 
dr r2 p dr 


describe the temperature variation in a convective element as it 

moves adiabatically (no radiative losses) through the surrounding 

layers. The adiabatic exponent I’, is given by 
32 — 246 — 3p? 

= 24—188 38? 

where f = p/p is the fraction of the total pressure contributed by 

gas pressure. We then set the actual temperature gradient to be, 

dT dTa  [min(—dTaa/dr, —dTya/dr) + dTraa/dr] 

dr dr +x") 

where x = F'/F con, max. Typically, we find that the pressure scale- 


height (which in a mixing length theory for convection is compara- 
ble to the distance convective elements travel) is a small fraction of 


T (43) 


» 44) 


3 Note that, in using the Schwarzschild rather than the Ledoux criterion for 
convective stability, we again ignore effects due to changes in the mean 
molecular weight with radius. As explained above, we believe that this is a 
good approximation. 
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the radius across the outer ~90 per cent of the quasi-star by radius. 
In the inner regions, where the pressure profile is quite flat, x is very 
small, and hence it is reasonable to assume that convection (even if 
it behaves differently from a mixing length theory) also succeeds in 
establishing the adiabatic gradient there. As in the analytic models, 
we note that the maximum flux that convection can transport is lim- 
ited by the condition that convective motions cannot be supersonic. 
Accordingly, we set 


Fon, max = Bes Pr (45) 


with the efficiency parameter 6 = 0.1. 

Equation (44) describes the limiting behaviours of the tempera- 
ture gradient in the different regions of a quasi-star. The adiabatic 
gradient applies if the region is dynamically unstable and the flux 
is less than the maximum convective flux; otherwise, the radiative 
gradient applies. These limiting cases accurately model the regions 
of high convective efficiency, in the denser and more opaque interior 
layers, and of high radiative efficiency, in the generally stable outer 
layer. In the transition region, however, our use of a flux limiter (the 
1 + x! term in equation 44) is only an approximation. A more ac- 
curate treatment would make use of mixing length theory. However, 
we have checked that different descriptions of the transition region 
do not affect appreciably the overall quasi-star structure. 


3.2 Integration method 


Equations (36), (37), (39) and (44), together with the definitions 
of the radiative and adiabatic temperature gradients, form a closed 
set of equations for the four unknowns (p, T, and mass), given 
a luminosity Lgy and stellar radius R,. To solve this system, we 
fix the photospheric temperature T,,, the black holes mass Mgy 
and the parameter œ, and guess the photospheric radius R, and the 
quasi-star mass M,. The luminosity is then Leg = 4nR20 Fis where 
o = ac/4. We determine the photospheric pressure and density us- 
ing the first equality in equation (18), together with the equation of 
state (equation 39). We then integrate inward from the known pho- 
tospheric conditions, adjusting our initial guesses until we match 
the two conditions: 


Lett = Lpu(10 x Reonai); (46) 


M(RBonai) =0. (47) 


We calculate Lgy (equation 3) with p and p evaluated at 10 x 
Rbgonai, Well beyond the inner region affected by the black hole grav- 
ity. For our numerical examples, œ is assumed to equal to 0.1, if 
not otherwise stated. We discuss the dependence of our results on 
a in Section 4.3. We explore two models for the opacity: the ‘toy’ 
opacity given by equation (15), which we used previously in the 
analytic models, and the full numerical opacity tabulated by MDOS5, 
which we refer to as ‘Pop III opacity’. Fig. 2 shows a comparison 
of these opacities for different densities of interest in the outermost 
radiative layers. 


3.3 ‘Toy’ opacity 


As a test, both of our analytic models and numerical scheme, we 
first compute structure models using the ‘toy’ opacity (equation 15). 
The toy opacity ignores density dependence, and in particular the 
contribution from bound-free and free—free absorption at T > 8 x 
10° K. As can be seen in Fig. 2, this is a reasonable approximation 
at low density, but is very poor at higher density. 


opacity [c mê gr’) 
a 


— "Toy" opacity 
al 


10 F -- -PoP Ill, p = 107"? gr cm™° 4 
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10° = 1 1 1 1 1 1 


44 46 48 5 5.2 5.4 
log10(T) [K] 

Figure 2. Opacity versus temperature. The Pop III opacity from MDO05 
is plotted as a function of temperature for three different values of matter 
density as shown in the legend. The range 107! < p < 1078 g cm~? spans 
the densities to be found in the outermost layers of quasi-stars. The solid line 
shows the analytic opacity given by equation (15). The figure shows how 
the Pop III opacity increases over the analytic fit as the density increases. The 
effect is especially evident around T = 10* K, where there is a bound-free 
peak due to hydrogen ionization. 
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Figure 3. Minimum photospheric temperature Tmin versus black hole mass 
in unit of solar masses mgp, for œ = 0.1. Numerical models with the ‘toy’ 
opacity (short-dashed line) are compared to results computed using the Pop 
IMI opacity (solid line). The analytic estimate, obtained by combining equa- 
tions (24) and (11), is shown as the long-dashed line. Static quasi-star models 
do not exist in the lower shaded region. 


As expected, we find that the analytic scalings are most reliable 
for high photospheric temperatures, and less reliable close to the 
temperature floor. We confirm that for a given (a, Mgu), there is 
indeed a minimum photospheric temperature, Tmin, above which 
the luminosity at the transition radius is sub-Eddington. Below Tmin, 
the whole radiative layer experiences a radiative force greater than 
the gravitational force and the quasi-star is bound to evaporate. In 
Fig. 3, the short-dashed line shows Tmin as a function of Mgy for 
numerical models computed using the toy opacity. The analytic 
estimate, obtained by combining equations (24) and (11), is shown as 
the long-dashed line. Numerically, we obtain a slightly higher Tmin- 
The discrepancy is higher for higher black hole masses. We note that 
for models computed using the toy opacity, the critical temperature 
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below which the entire radiative layer becomes super-Eddington 
(our definition of Tmin) is close to the temperature at which any 
point in the radiative layer becomes super-Eddington. This simple 
behaviour, which occurs because the toy opacity is constant in most 
of the radiative layer before dropping off steeply and monotonically 
near the photosphere, does not carry over to models computed using 
full Pop III opacity. 


3.4 Pop III quasi-stars 


To extend these results, we compute numerical quasi-star struc- 
tures with realistic opacities. We assume that quasi-stars have a pri- 
mordial (metal-free) nuclear composition, and use the opacity table 
calculated by MDOS. This table covers the density range 107!6 < 
p(gcm~*) < 107? for temperatures 63 < T [K] < 4 x 10*. To com- 
pute the interior structure of the quasi-star — where temperatures 
attain values much higher than 4 x 10+ K — we analytically extend 
the opacity table assuming that the excess opacity above the electron 
scattering value has a Kramers form, i.e., 


K(p, T) = CPT + Kes. (48) 


We fix the function C(p) to match smoothly on to the tabulated 
opacity at the highest temperature point. As expected, Thomson 
scattering is the only significant source of opacity deep within the 
interior. 

As is evident from our Fig. 2, or from fig. 4(a) of MDOS, the 
opacity deviates substantially from our toy model around T ~ 104 
K for p > 107? gcm~?. The prominent enhancement of the opacity 
with increasing density is due to neutral hydrogen absorption, H + 
hv > H+ + e. 

The differences between the toy and full Pop II opacities are 
significant only at low temperatures. In the deep interior of the quasi- 
star, the opacity is constant, as it was for the analytic opacity, and we 
again find a convective zone that is described by an n = 3 polytrope 
for which pg/p, = 7.1m;,°°. 

In the outer region, the flux is transported by radiative diffusion. 
The radiative layer covers between ~50 and ~10 per cent of the 
quasi-star by radius, becoming thinner as the photospheric temper- 
ature decreases towards Tmin (for a fixed mpy) or as the black hole 
becomes more massive (for a fixed Tpn). In contrast, the per cent of 
mass in the outer radiative layer increases as the limiting temper- 
ature is approached. For mgy = 300, it goes from ~1 per cent at 
Tpn = 10* to ~18 per cent around Tmin, While for mgy = 7 x 10° it 
shows a similar behaviour but never exceeds 2 per cent. 

The dramatic differences between structures computed using the 
real opacity, and those based on the toy opacity, are almost exclu- 
sively confined to the radiative layer. In Figs 4 and 5, we show 
examples of the behaviour of the opacity and the local Eddington 
ratio, Lraa/Leaa With Leda = 4GM(r) c/x(r), for models computed 
with four different photospheric temperatures and constant black 
hole mass. The plots focus on the radiative zone. For clarity, we 
plot temperature on the x-axis, since it spans more than an order of 
magnitude while the radius increases only by tenths of a per cent. As 
Tpn decreases, the opacity increases everywhere and the growth of 
the bound-free peak around T = 10* K becomes particularly promi- 
nent. Correspondingly, a peak in the Lyaq/Leaa ratio forms. This peak 
first becomes super-Eddington for T ~ 10* K, and steadily grows 
and expands as the photospheric temperature drops until the entire 
radiative zone is super-Eddington at Ty, = Tmin ~ 4 x 10° K. Note 
that also for T,, ~ 10* K, a narrow super-Eddington region exists 
around T ~ 3 x 10* K. This is also true for higher photospheric 
temperatures. We comment on this point at the end of this section. 
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Figure 4. The opacity as a function of temperature in the outer layers of 
four quasi-star models with Mgy = 300 Mọ and « = 0.1, computed using 
Pop III opacities. The photospheric temperatures Tph of the models are shown 
in the legends. On the right-hand panel, we adopt a logarithmic scale for clar- 
ity, while the plot on the left-hand panel uses a linear scale. The temperature 
at the transition between the convective and radiative zones is denoted by the 
corresponding vertical lines. Note the rise in the bound—free peak at around 
T ~ 10* K (marked on the left-hand panel by an arrow), which becomes 
dominant as the photospheric temperature drops and the photospheric den- 
sity increases (see also Figs 7 and 8). Note also the steep drop in opacity at 
the photospheres. 
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Figure 5. Lraq/Leaa — 1 for the four models shown in Fig. 4. The solid 
horizontal line marks the border between the super-Eddington (above) and 
the sub-Eddington (below) regions. Note that for Tph ~ 10* K, the super- 
Eddington zone is very narrow and confined around T ~ 3 x 10* K. For 
Tph < 10* K, a peak appears around T ~ 104 K (marked by an arrow in the 
left-hand panel) and expands, eventually encompassing the whole radiative 
layer (dotted line, right-hand panel). 


Fig. 4 clearly shows the steep drop in opacity at the photosphere, 
and how it becomes more vertiginous as the quasi-star cools down 
to Tmin. This feature ensures (as for the ‘toy’ opacity) the presence 
of a minimum temperature, below which no hydrostatic solutions 
can be found (Fig. 3). As Tpn decreases, the boundary condition at 
the photosphere (equation 18) implies that the photospheric pres- 
sure passes through a minimum and then begins to increase. The 
photospheric radius and the transition radius Ry, then stop increas- 
ing and start sinking rapidly inward in order to maintain hydrostatic 
balance. Eventually, no hydrostatic solution is possible. 

Despite the complicated behaviour of the opacity, the ana- 
lytic power-law scalings between quasi-star mass, photospheric 
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Figure 6. The quasi-star mass is plotted as a function of the photospheric 
temperature (solid line) for models with mgy = 300 Mọ and a = 0.1, 
computed with full Pop III opacities. The dotted line shows the analytic 
scaling m, « To , Obtained assuming that £w = K = 1. The dashed 
line is also derived from the analytic relation given in equation (11), but 
with numerical values for £y and K used instead of assuming that these 
factors are exactly unity. The star shows the mass of a (unphysical) numerical 
model computed with no radiative zone. This closely matches the analytic 
prediction, and indicates that the offset at high temperatures is due to the 
presence of a radiative zone. 


temperature and black hole mass (equation 11) remain reliable at 
sufficiently high temperatures, although they suffer an offset in nor- 
malization correlated with the radial thickness of the radiative layer. 
To illustrate this, Fig. 6 shows the numerically computed m, versus 
Tpn for models with Mgy = 300 Mo. The analytic scaling, 


m, x Ta (49) 


holds well down to Tp, ~ 7 x 10° K, though there is an offset in 
the normalization. This offset is due to the presence of the radiative 
layer. A fully convective model at these high photospheric temper- 
atures would lie on the analytic track: an example of such a model 
calculated for T,, = 9000 is plotted in Fig. 6 as a star. For higher 
mpu, the offset is reduced as the radiative layer gets thinner. At 
lower temperatures, m, deviates slowly from the power law until 
Tmin is approached, whereupon m, drops more steeply. This is the 
catastrophic shrinking of the quasi-star prior to evaporation. In this 
regime, the analytic scalings do not hold. 

Although the integrated properties of the numerical Pop III opac- 
ity models resemble those derived analytically and numerically us- 
ing the toy opacity, the presence of locally super-Eddington fluxes 
for photospheric temperatures T,, > Tmin has substantial implica- 
tions for the derived structures. Figs 7 and 8 show the radial profiles 
of temperature and density for the models shown in Figs 4 and 5. In 
hydrostatic equilibrium, the combination of equations (36) and (39) 
yields an expression for the density gradient, 


kT do GM(r) dp pk aT 
= p 


(50) 


ud r2 dr u dr` 


In regions where the radiative force — dp,/dr substantially exceeds 
the gravitational force, the density gradient becomes positive and a 
local density inversion forms. At the same time, the temperature gra- 
dient steepens. In turn, the increase in the density and the decrease 
in temperature enhance the contribution from bound—free hydro- 
gen absorption, raising the opacity further above Kes. This general 


Figure 7. Radial profiles of the interior temperature for the four models 
plotted in Fig. 4. For each model, the transition radius Ry is marked by the 
corresponding vertical line. 
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Figure 8. The same as Fig. 7, but for density. As Tpn decreases, the radiative 
layer becomes denser. 


behaviour becomes more marked as the photospheric temperature 
drops, and the radiative layer becomes denser and cooler. 

The existence of density inversions in our models immediately 
raises the question of whether such structures are stable. There are 
two conceptually distinct concerns. In one dimension, the hydro- 
static configurations we have computed are self-consistent, provided 
that the super-Eddington regions lie beneath the photosphere. An 
inward-directed force due to the gas-pressure gradient compensates 
for the imbalance between radiation pressure and gravity. However, 
there could exist a second solution in which the super-Eddington 
flux drives a wind from the quasi-star. If such solutions exist, the 
hydrostatic solution could be unstable even in one dimension, lead- 
ing to mass-loss at temperatures well above Tmin. In two dimensions, 
itis even more doubtful that a density inversion would be stable. The 
resulting instabilities could lead to lateral density contrasts and non- 
magnetized photon bubbles of the type analysed by Shaviv (2001). 

We have not had notable success in analytically estimating the 
mass-loss rate that might occur due to locally super-Eddington fluxes 
in the quasi-star, and our numerical scheme is not suited to tackle 
the problem. This failure is hardly surprising, since the problem of 
what happens when a star develops a limited super-Eddington zone 
has been studied extensively in the context of luminous blue vari- 
ables such as n Carina (for a review, see Humphreys & Davidson 
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1994), apparently without any definitive theoretical resolution be- 
ing attained. One possibility is that a large-scale circulation pat- 
tern develops, superficially resembling convection, and the star suf- 
fers no mass-loss at all (Owocki, Gayley & Shaviv 2004). Another 
possibility, empirically favoured for LBVs, is that episodic mass- 
loss occurs. 

For what follows, we conjecture that since the super-Eddington 
zones in our models arise due to high densities (and vanish if the 
density is reduced), quasi-stars with Pop III opacities may develop 
outer regions in which there is circulation but little or no mass-loss. 
As the width of the super-Eddington region becomes comparable 
to the radiative layer thickness, however, the extended acceleration 
zone is likely to permit a strong wind to develop. We therefore 
adopt Tmin aS an estimate for the minimum temperature a quasi-star 
can sustain before evaporating. That there are uncertainties in this 
identification should be very obvious. 


4 CO-EVOLUTION OF BLACK HOLES 
AND QUASI-STARS 


The properties of a quasi-star will change as the black hole grows 
and, possibly, as the quasi-star itself accretes matter from its environ- 
ment. The thermal time-scale in a quasi-star interior is sufficiently 
short that the structure can adjust quasi-statically. Therefore, we 
can model the co-evolution of the black hole and envelope as a se- 
ries of equilibrium models, as long as the photospheric temperature 
exceeds the floor associated with the opacity crisis, Ty, > Tmin- 

We will consider two models for the quasi-star’s interaction 
with its environment. In the first, we will assume that the quasi- 
star accretes matter at a constant rate, which we parametrize as 
0.1719.1 MO yr—!. This growth rate is consistent with the BVR sce- 
nario for black hole growth in pre-galactic haloes. In the second, 
we will assume that the quasi-star has a fixed mass, which we 
parametrize in units of 10° Mọ. We normalize time in years and 
assume, again motivated by the BVR arguments, that the initial 
seed black hole mass is much smaller than the final mass attained 
prior to dissolution of the quasi-star. 


4.1 Before the opacity crisis 


As long as Tp, > Tmin, the growth rate of the black hole is set by 
the Eddington limit of the quasi-star, implying 


fitga = 2.5 x 10-8e5 1m, Mo yr, (S1) 


where € = 0.1é€o, is the accretion efficiency and we have taken 
fy = K = 1. The accretion efficiency relates the luminous output 
due to accretion with the rate of growth of the black hole mass. 
This efficiency factor is distinct from the Bondi efficiency factor a 
defined in equation (3), and could be smaller than the standard value 
of 0.1 if energy is lost mechanically in a jet that does not couple to 
the envelope. When the quasi-star grows at a steady rate, we have 
m, = 0.179; ty, and the black hole mass grows according to 


mpy = 1.2 x 10 eG {rit ity, = 1.2 x 107 eji mgin. (52) 


If the quasi-star has a fixed mass, then the black hole grows linearly 
with time: 


mpu = 2.5 x 10779 i msstyr. (53) 


Equations (11) and (52) imply that the photospheric temperature de- 
creases as the black hole grows. The quasi-star thus evolves towards 
the opacity crisis, reaching it when Tph = Tmin = 4000 Tm,4 K. For 
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the case of steady quasi-star growth, this occurs when the quasi-star 
mass reaches 


m, = 1.8 x 10° oP mast aa 
and the black hole mass reaches 
mpy, = 3.9 x 10°) mn 0, OTA, ies 


m 


For a fixed envelope mass, the opacity crisis is reached when 
mpno = 1.9 x 10*09 mE T 4”. (56) 
4.2 Numerical results for quasi-star evolution 


To test the above results, we numerically computed explicit evolu- 
tionary sequences by solving the equations, 


dMgy Legy 
ay BE. 57 
dt ec? On 
oie = M (58) 
dt si x9 


with Ley(M., Mgu, œ) calculated from the model with no assump- 
tion as to the luminosity relative to the Eddington limit. The analytic 
evolution tracks for quasi-stars depend only on the assumption that 
the luminosity of the quasi-star is close to the Eddington limit ap- 
propriate to the total mass. Our numerical results verified that this 
is a very good approximation. Accordingly, to estimate the maxi- 
mum mass a black hole can grow to within a quasi-star it suffices to 
use the analytic growth track given by equation (52). The analytic 
estimate of the minimum temperature is less accurate. We therefore 
combine the analytic growth track with the numerically computed 
minimum photospheric temperature to derive the final mass. 

In Fig. 9, we plot analytic tracks corresponding to envelope accre- 
tion rates of 0.1 Mọ yr! and 1 Mo yr!. As the black hole grows, 
the tracks move towards the right-hand side of the plot, eventually 
crossing into the forbidden region set by the minimum photospheric 
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Figure 9. Envelope mass versus black hole mass for static solutions at the 
minimum photospheric temperature. The dashed line is for the ‘toy’ opacity 
model, the solid lines are for the numerical opacity. Static solutions are ex- 
cluded in the lower shaded regions. The lighter shaded region is computed 
assuming that a = 0.1, the darker shaded region is for œ = 0.05. Super- 
imposed on the figure are the evolutionary tracks (equation 52). The upper 
track is for an accretion rate on to the envelope of M,= 1Mọ yr™!. The 


lower track is for M, = 0.1 Mo yr, 
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temperature. The figure shows that the final black hole mass is higher 
for higher accretion rates on to the envelope, because the effective 
temperature decreases more slowly, delaying the final dissolution 
at Tmin. We find that quasi-stars are indeed an efficient channel for 
‘growing’ intermediate-mass black holes. For a = 0.1 and mo1 > 1, 
the final black hole mass is predicted to be at least a few thousand 
solar masses. 


4.3 Scaling of results with a 


The parameter œ, which accounts for inefficiencies in the accretion 
flow within the Bondi radius, is largely unknown. The numerical 
results we have calculated are mostly for a = 0.1. Analytically, 
there is a simple scaling with a, valid for m, >> mgu. In this limit, 
the only dependence of the quasi-star structure on the black hole 
mass enters via equation (3) for the luminosity, Leu « M?,0. Any 
solution in which the combination M = Mpya!/? is constant should 
then have the same m, and T,,, along with the same radial structure 
well outside the Bondi radius. If œ were to be smaller than our 
assumed value of 0.1, this would allow larger black holes to grow 
within the same envelopes. 

Numerically, the expectation that the results depend only on M 
is borne out for mgy > 10%, irrespective of the photospheric temper- 
ature. For mgy < 10°, it is valid only for photospheric temperatures 
well in excess of the minimum temperature. For low black hole 
masses and photospheric temperatures approaching the minimum 
value, the black hole mass cannot be ignored in the equation of hy- 
drostatic equilibrium, and the resulting structure depends separately 
on mpy and æ. In Fig. 9, we plot numerical results for the minimum 
photospheric temperature computed with a = 0.1 and a = 0.05. 
As is evident from the figure, the offset between these curves is not 
constant. However, for the higher masses that are of interest when 
determining the maximum black hole mass that can be attained, the 
a!/? scaling is quite accurate. 


4.4 Post-opacity crisis: dissolution of the quasi-star 


Our analysis of the junction between the radiative layer and convec- 
tive zone indicates that no static solutions exist for Ty, < Tmin- In 
contrast to a red giant, protostar or Thorne—Zytkow object, where 
feedback allows the energy source or envelope to adjust so that the 
photosphere follows the Hayashi track, no stable feedback appears 
to exist here. Once a quasi-star reaches Tin, the plummeting opacity 
causes the convective zone to release radiation at a super-Eddington 
rate. The deflation of the convective zone increases the rate of ac- 
cretion on to the black hole, leading to a runaway. The only ways 
to avoid the destruction of the quasi-star would be to decrease the 
mass of the black hole or to increase the mass of the envelope at an 
unrealistically high rate. 

A detailed analysis of the mass-loss process is beyond the scope 
of this paper, but it is easy to show that the black hole is unlikely 
to grow much while the quasi-star is evaporating. Once mass-loss 
starts in earnest, the wind will carry away nearly all the energy 
released by the black hole. Assuming that the wind speed is of 
order of the escape speed from the quasi-star, vw ~ (GM, /R,)'/ 2. 
this implies that the quasi-star evaporates when the energy liberated 
by accretion equals the binding energy. (Although the quasi-star is 
mainly radiation-pressure dominated, we assume that the binding 
energy is enhanced by rotation.) The mass accumulated by the black 
hole during the dissolution phase is then 

2 
AMpu ~ u Wes hga Uy Tra MO, (59) 
*0 


which is negligibly small compared to mpgy,o. Therefore, we may re- 
gard mpguy,o as the maximum mass attainable by a black hole growing 
inside a quasi-star envelope. 


5 DISCUSSION AND CONCLUSIONS 


In this paper, we have studied the structure and evolution of quasi- 
stars, rapidly accreting black holes embedded within massive gas 
envelopes. We find that, for any black hole mass, there is a minimum 
photospheric temperature below which rapid dissolution of the en- 
velope is inevitable. Both analytic and numerical models, computed 
using Pop III opacities, suggest that this minimum temperature is 
around 4000-5000 K. If quasi-stars are implicated in the forma- 
tion of seeds for supermassive black holes in pre-galactic haloes, 
as suggested by BVR, this floor temperature implies that the most 
luminous quasi-stars would emit most strongly in the rest-frame 
near-infrared. At typical redshifts of z ~ 10, the observed spec- 
trum would peak at à ~ 10 um. We defer discussion of the possible 
cosmological density of such sources, and hence their observability 
with future facilities such as the James Webb Space Telescope, to a 
subsequent paper. 

We have also studied the evolution of quasi-stars in simple sce- 
narios for the mass growth of their envelopes. Generically, as the 
black hole mass grows the photospheric temperature falls, un- 
til eventually the limit imposed by the behaviour of the Pop III 
opacity is reached. Unveiling of the black hole appears inevitable 
long before it succeeds in accreting much of the envelope. Both 
analytic and numerical estimates suggest that seed black holes 
with masses between a few 10° and 10+ Mo are plausible out- 
comes of this scenario. The efficiency of the black hole accretion 
process that powers quasi-stars (parametrized here by œ) is un- 
known — if the efficiency is low (due, e.g. to polar outflows that 
couple poorly to the envelope) then larger seed black holes are 
possible. 

Although our numerical integrations of quasi-star structure con- 
firm many aspects of the analytic model, they also reveal complex 
behaviour in the radiative zone immediately beneath the photo- 
sphere. Narrow regions in which the flux is locally super-Eddington 
develop for temperatures significantly in excess of the minimum 
temperature at which the entire radiative zone becomes super- 
Eddington. In our hydrostatic models, these zones are character- 
ized by density inversions. Interpreting this structure is tricky, 
since such density inversions may well be subject to one- or two- 
dimensional instabilities whose ultimate resolution is unknown. 
Hydrodynamic studies will be needed to determine the extent of 
mass-loss that may occur at temperatures above the theoretical floor 
value. 

Finally, we note that here we have focused on black holes em- 
bedded within truly primordial gas at high redshift. Small amounts 
of metal pollution would act to increase the opacity in the radiative 
zone, altering the minimum temperature and possibly increasing 
the likelihood of mass-loss due to the formation of dust. Changes to 
the interior structure would be smaller. In particular, the density at 
the base of the radiative zone is so low (p ~ 107° g cm~?) that elec- 
tron scattering will continue to dominate the opacity for T > 10° K, 
even in the presence of pollution. Of greater concern is the fact that 
metal-enriched gas in the halo outside the quasi-star would be more 
susceptible to fragmentation and star formation. Too much star for- 
mation would reduce the rate of mass accretion on to the quasi-star 
below the values (~0.1 Mo yr”! or higher) that we have assumed, 
resulting in smaller final black hole masses. It is also possible that 
generically similar structures could form at lower redshift whenever 
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small black holes encounter very high rates of gas inflow. Structures 
similar to those we have described could allow stellar remnants at 
the centre of merging galaxies to grow significantly via a transient 
quasi-star stage. 
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